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Abstract. In this paper we study the spectrum S of the infinite Feinberg- 
Zee random hopping matrix, a tridiagonal matrix with zeros on the main 
diagonal and random il’s on the first sub- and super-diagonals; the study of 
this non-selfadjoint random matrix was initiated in Feinberg and Zee {Phys. 
Rev. E 59 (1999), 6433-6443). Recently Hagger {Random Matrices: Theory 
AppL, 4 1550016 (2015)) has shown that the so-called periodic part S,r of S, 
conjectured to be the whole of S and known to include the unit disk, satisfies 
p“^(E.,r) C Titt for an infinite class S of monic polynomials p. In this paper 
we make very explicit the membership of S, in particular showing that it 
includes Pm{\) = At/m-i(A/2), for m > 2, where U„{x) is the Chebychev 
polynomial of the second kind of degree n. We also explore implications of 
these inverse polynomial mappings, for example showing that is the closure 
of its interior, and contains the filled Julia sets of infinitely many p £ <S, 
including those of Pm, this partially answering a conjecture of the second 
author. 


Mathematics Subject Classification (2000). Primary 47B80; Secondary 37F10, 
47A10, 47B36, 65F15. 


Keywords, random operator, Jacobi operator, non-selfadjoint operator, spec¬ 
tral theory, fractal, Julia set. 



2 


Simon Chandler-Wilde and RafFael Hagger 


1. Introduction 

In this paper we study spectral properties of infinite matrices of the form 

(■■ \ 


Ar = 


0 1 
CO 0 1 

Cl 0 


V 


( 1 . 1 ) 
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where c G := {±1}^ is an infinite sequence of ±l’s, and the box marks the 
entry at (0,0). Let denote the linear space of those complex-valued sequences 
(/) : Z —>• C for which \\4>\\2 '■= < oo, a Hilbert space equipped 

with the norm |j • II 2 . Then to each matrix A^. with c G corresponds a bounded 
linear mapping £‘^ ^ (A^ which we denote again by Ac, given by the rule 

(Ac</^)m ~ ^7n4^m—l “t” ^m+l, ^ ^ 

for (j) G 1“^. 

Following we will term (1.1) a Feinberg-Zee hopping matrix. Further, in 
the case where each Cm is an independent realisation of a random variable with 
probability measure whose support is { — 1,1}, we will term Ac a Feinberg-Zee 
random hopping matrix, this particular non-selfadjoint random matrix studied 
previously in [liiiniiHiiiiiEiigiiiiisiiieiini.Q The spectrum of a realisation Ac 
of this random hopping matrix is given, almost surely, by (e.g., 0) 


spec He = E := (J specHf,. 
beci 


( 1 . 2 ) 


Here specHf, denotes the spectrum of Af, as an operator on Note that (1.2) 
implies that E is closed. 

Equation (O holds whenever c G 11 is pseudo-ergodic, which means simply 
that every finite sequence of ±l’s appears as a consecutive sequence somewhere in 
the infinite vector c; it is easy to see that c is pseudo-ergodic almost surely if c is 
random. The concept of pseudo-ergodicity dates back to [5], as do the arguments 
that (1.2) holds, or see [1] for <[H derived as a special case of more general limit 
operator results. 

Many of the above cited papers are concerned primarily with computing up¬ 
per and lower bounds on E. A standard upper bound for the spectrum is provided 
by the numerical range. It is shown in 0 that, if c S H is pseudo-ergodic, its nu¬ 
merical range W{Ac), defined by W{Ac) '■= {{Ac(t>,4>) '■ (t> G £'^, \\(t >\\2 = 1}, where 


^ These random hopping matrices appear to have been studied initially in m> in which paper 
the first superdiagonal is also a sequence of random ±l’s. But it is no loss of generality to restrict 
attention to matrices of the form as the case where the superdiagonal is also random can 

be reduced to 113 by a simple gauge transformation; see m or [ 4 ] Lemma 3.2, Theorem 5.1]. 
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(•, •) is the inner product on is given by 

W{Ac) = A:={x + iy: x,y G K, |a;| + \y\ < 2}. 


(1.3) 


This gives the upper bound that S C A, the closure of A. Other, sharper upper 
bounds on S are discussed in Section [2] below. 

This current paper is related to the problem of computing lower bounds for E 
via If & G n is constant then At, is a Laurent matrix and specH{, = [—2, 2] if 
bm = 1, while specHf, = i[—2,2] if = —1; thus, by ( |1.2[ ), tti := [—2, 2]Ui[—2, 2] C 
E. Generalising this, if 6 G O is periodic with period n then specHf, is the union of 
a finite number of analytic arcs which can be computed by calculating eigenvalues 
of n X n matrices (see Lemma 2.2 below). And, by (1.2), 7r„ C E, where 7r„ is the 


union of specA^ over all b with period n. This implies, since E is closed, that 


^TT ■- A E, 


(1.4) 


where := U„gN7r„. 

We will call E,r the periodic part of E, noting that [3] conjectures that equality 
holds in (1.4), i.e. that tToo is dense in E and E,r = S. Whether or not this holds is 
an open problem, but it has been shown in [S] that tToo is dense in the open unit 
disk D := {A G C : |A| < 1}, so that 


C E^ C E. 


(1.5) 


For a polynomial p and S' C C, we dehne, as usual, p{S) := {p{X) : A G S} 
and p~^{S) := {A G C : p{X) G S}. (We will use throughout that if S is open then 
p~^{S) is open (p is continuous) and, if p is non-constant, then p{S) is also open. 


e-g-, [m Theorem 10.32].) The proof of (1.5) in depends on the result, in the 
case p{X) = A^, that 


P ^(tToo) C tToo, so that also p ^(E,n.) C E,n 


( 1 . 6 ) 


This implies that S„ C tToo, for n = 0,1,..., where So := [—2,2] and Sn := 
p“^(S„_i), for n G N. Thus U„ gPjSn, which is dense in D, is also in tToq, giving 


(1.5). 


Hagger m makes a large generalisation of the results of [S], showing the 
existence of an infinite family, S, containing monic polynomials of arbitrarily high 


degree, for which (1.6) holds. For each of these polynomials p let 


u{p) ■- U P~ 


(1.7) 


n—1 


(Here p ^(S) := p ^{p ^(S)), p ^(S) := p ^{p ^(S)), etc.) Hagger [T7] observes 
that, as a consequence of (1.5) and (1.6), U{p) C E,r- He also notes that standard 


results of complex dynamics (e.g., El Corollary 14.8]) imply that J{p) C U{p), 
so that J{p) C E,r; here J{p) denotes the Julia set of the polynomial p. (Where 
p^(A) := p{p{X)), p^(X) := p(p^(A)), etc., we recall [TT] that the filled Julia setK{p) 
of a polynomial p of degree > 2 is the compact set of those A G C for which the 
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sequence (p"'(A))„gN, the orbit of A, is bounded. Further, the boundary of K{p), 
J{p) := dK{p) C K{p), is the Julia set of p.) 

The definition of the set S in while constructive, is rather indirect. The 
first contribution of this paper (Section is to make explicit the membership of 
S. As a consequence we show, in particular, that € 5, for m = 2,3 ,..., where 
Pm (A) := AC/m -1 (A/2), and Un is the Chebychev polynomial of the second kind 
of degree n [T]. 

The second contribution of this paper (Section]^ is to say more about the 
interior points of Previous calculations of large subsets of tToo, precisely calcu¬ 
lations of TTji for n as large as 30 laii], suggest that Stt fills most of the square A, 
but int(I]^), the interior of is known only to contain D. Using that the whole 
family {Pm : m > 2} C 5, we prove that (—2,2) C int(I] 7 r). This result is then 
used to show that is the closure of its interior. Using that p“^(]D)) C E^, for 
p G S, we also, in Section 3.2 construct new explicit subsets of E^ and its interior; 
in particular, extending (1.5), we show that aD C E.,r for a = 1.1. 


In the final Section of the paper we address a conjecture of Hagger m 
that, not only for every p £ S is J{p) C U{p) (which implies J{p) C E^), but also 
the filled Julia set K{p) C U{p). This is a stronger result as, while the compact set 
J{p) has empty interior [TTl Summary 14.12], K{p) contains, in addition to J{p), 
all the bounded components of the open Fatou set F{p) := C \ J[p). We show, 
by a counter-example, that this conjecture is false. But, positively, we conjecture 
that K{p) C Ett for all p £ S, and we prove that this is true for a large subset of 
S, in particular that K{Pm) C E^ for m>2. 

The results in this paper provide new information on the almost sure spec¬ 
trum E D E,n. of the bi-infinite Feinberg-Zee random hopping matrix. They are 
also relevant to the study of the spectra of the corresponding finite matrices. For 
n S N let Vn denote the set of n x n matrices of the form (1.1), so that Vi := {(0)} 
and, for n > 2, U„ := : k = {ki, ..., fc„) e {il}"}, where 

/ 0 1 \ 


aJ”) := 


0 1 
fci 0 


V 


( 1 . 8 ) 


kn-l 0 / 

(This notation will be convenient, but note that A["^ is independent of the last 
component of k.) Then spec A^"^ is the set of eigenvalues of the matrix A[”^ Let 


OO 

cr„ := spec A, for n S N, and a^o ■= [J cr„, (1.9) 

AgVti n=l 

so that (Too is the union of all eigenvalues of finite matrices of the form ( |1.8| ). 
Then, connecting spectra of finite and infinite matrices, it has been shown in [4] 
that cr„ C TT2n+2, for n € N, so that a^o C C E^r. Further, [T^ shows that is 
dense in tToo, so that E^^. = In Section [XT] we build on and extend these results, 
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making a surprising connection between the eigenvalues of the finite matrices (1.8) 
and the spectra of the periodic operators associated to the polynomials in S. The 
result we prove (Theorem |3.8[), is key to the later arguments in Section]^ 


2. Preliminaries and previous work 

Notions of set convergence. We will say something below about set sequences, 
sequences approximating S and from above and below, respectively. We will 
measure convergence in the standard Haussdorf metric d(-,-) [TBl Section 28] (or 
see El) on the space C*" of compact subsets of C. We will write, for a sequence 
(Sn) C and S € C^, that Sn \ S if d{Sn, 5) —>■ 0 as n —)■ oo and S C Sn for 
each n; equivalently, Sn \ S as n ^ oo if S C Sn for each n and, for every e > 0, 
Sn C S + eD, for all sufficiently large n. Similarly, we will write that Sn S if 
d{Sn, S') —?► 0 as n —>■ oo and S„ C S for each n; equivalently, S„ S if S„ C S 
for each n and, for every e > 0, S C S„ + eD, for all sufficiently large n. The 
following observation, which follows immediately from m Proposition 3.6] (or see 
[S] Lemma 4.15]), will be useful: 

Lemma 2.1. If Si C S 2 C ... C C are closed and Soo ■= is bounded, then 

Sn Soo, as n —>■ 00. 


Spectra of periodic operators. We will need explicit formulae for the spectra of 
operators Ac with c G in the case when c is periodic. For k = (fci,...,fc„) € 
{±1}”, let denote Ac in the case that Cm+n = c„, for m S Z, and Cm = ^m, 
for m = 1,2, ...,n. For n G N let denote the order n identity matrix, 
the n X n matrix which is zero except for the entry 1 in row n, column 1, and 
let Rn denote the transpose of i?„. For n e N, fc G {±1}", and ip G R, let 
ak{ip) := A^^'^ + e~^‘^Rn + knO^'^R^- The following characterisation of the spectra 
of periodic operators is well-known (see Lemma 1 and the discussion in |16jL 

Lemma 2.2. For n G N and k G {±1}", 

specH^®'^ = {A G C : det(afe((^) — A/„) = 0 for some (p G [0, 27r)}. 


Key to our arguments will be an explicit expansion for the determinant in 
the above lemma, expressed in terms of the following notation. For n G N, fc = 
(ki ,..., kn) G {±1}", and A G C, let 


9fc(A):= 


A 1 

ki 

1 

kn A 


( 2 . 1 ) 
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For j S Z and A G C, let k[i : j) := [ki,kj), for 1 < i < j < n, and define 

A, ifi-j = l, 


Qk{i:j) ('^) • 


1 , 

0 , 


if J - j = 2, 
if * - J = 3. 


( 2 . 2 ) 


Then, for n G N and k G {±1}", expanding the determinant (2.1) by Laplace’s 
rule by the first row and by the last row, we see that 

^fc(A) = Af/fc^2:n)(A) ^l^fc(3:n)(A) = A^7/i;(l:n—1) (A) ^nQ/c(l:n—2) (A). (2-3) 

The following lemma follows easily by induction on n, using (2.3). The bounds on 
qk stated are used later in Corollary |5.2| 


Lemma 2.3. If k = (fci,...fc„) G {±1}”, for some n G N, then qk is a monic 
polynomial of degree n + 1, and qk is even and qk(0) = ±1 ifn is odd, qk is odd and 
<?fe(0) = 0 ifn is even. Further, |gfc(A)| > |(?fc(i:„_i)(A)| + l, |qfc(A)| > |gfe( 2 :„)(A)| + l, 
and |(7fc(A)| > n + 2, for |A| > 2. 

For n G N let J„ denote the n x n flip matrix, that is the n x n matrix with 
entry Si^n+i-j in row i, column j, where is the Kronecker delta. Then 
so that (det J„)^ = 1. For k = (fci,...,fc„) G {il}”, let k' := kJn = (fc„,...,fci). 
The first part of the following lemma is essentially a particular instance of a general 
property of determinants. 

Lemma 2.4. If k G {±1}", for some n G N, and I = fc', then qk = qi; if (■ = —k, 
then 

qi{X) =i~’^~^qkiiX). (2.4) 


Proof Suppose first that I = k!. Then qk{X) given by (2.3) is the determinant 
of a matrix A, and q£(X) is the determinant of J„+iA^ J„+i, so that g£(X) = 
(det J„+i)^ det Al = qk{X). That (2.4) holds if ^ = —k can be shown by an easy 

□ 


induction on n, using ( |2.3[ ), or directly by a gauge transformation. 

Flere is the announced explicit expression for the determinant in Lemma |2.2[ 
Lemma 2.5. [171 Lemma 3] For n gN, k G {±1}", A G C, and (p G [0, 27r), 

n 

(-1)” det(afc(</ 5 ) - Xln) = Pk{X) - e“^ kj - e”'"^, 

1=1 

where pk is a monic polynomial of degree n given by 

( f) Pk{X) = qk(l:n — l){ A) knqk{2:n—2){ A). (2-^) 

Further, pk is odd (even) ifn is odd (even). 

Since, from the above lemmas, pk is odd (even) and qk even (odd) if n is odd 
(even), (2.5) implies that 

PkiX') = qk{l:n—l){X) knqk{2:n—2){X') (^-S) 

= <Zfe(l:n-l)('^) + 9fc(2:n)(-^) ~ '^9fe(2:n-l)('^)j (2-7) 
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this last equation obtained using (2.3). The following lemma, proved using these 
representations, makes clear that many different vectors k correspond to the same 
polynomial pk- 


Lemma 2.6. If k = {ki, ...kn) G {±1}", for some n G N, and £ = k' or £ is a cyclic 
permutation of k, then pk = pe- If £ = —k then pg{X) = i“"pfc(iA). 


Proof. Using ( |2.6[ ) and ( |2.3| ) we see that 

PiiX) Pk(.X) = ^£(l:n—T ^nQk{2:n—2) ('^) InQi(2:n—2)(.X) 

= {Qe(2:n-l){X) — qk(l:n-2)iX)) — ^l9^(3:n-l) (^) 

+ kn-iqk{l:n-3)iX) + ^n'Zfc(2:ri-2) (A) — £nqi{2-.n-2){X) ■ 

If £ is a cyclic shift of A:, i.e., £j = %-i, j = 2, ...,n, and £i = kn, then the 
right-hand side is identically zero. Thus pg = if £ is a cyclic permutation of k. 

If £ = fc' then that pk = pg follows from (2.7) and Lemma 2.4 If £ = —k then 
that pg{X) = i“"pfc(iA) follows from (2.6) and Lemma [2.4| □ 


Call k G {±1}" even if ]Xj=i = 1) and odd if 117=1 


Corollary 5], it is immediate from Lemmas 2.2 and 2.5 that 


U=i 


-1. Then [m 


specH^®"^ = 2 , 2 ]), if k is even, specH^®'^ = p^ ^(i[— 2 , 2 ]), if k is odd. 

( 2 . 8 ) 

Complex dynamics. In Section below we show that filled Julia sets, K{p), of 
particular polynomials p, are contained in the periodic part of the almost sure 
spectrum of the Feinberg-Zee random hopping matrix. To articulate and prove 
these results we will need terminology and results from complex dynamics. 

Throughout this section p denotes a polynomial of degree > 2. We have 
defined above the compact set that is the filled Julia set K{p), the Julia set J(p) = 
dK{p) C K{p), the Fatou set F{p) (the open set that is the complement of J{p)), 
and the orbit of z G C. It is easy to see that, if z ^ K{p), i.e., the orbit of z is not 
bounded, then p"(z) —>■ oo as n —)■ oo, i.e., z G Hp(c»), the basin of attraction of 
infinity. We call S C C invariant if p{S) = S, and completely invariant if both S 
and its complement are invariant, which holds iff p~^{S) = S. Clearly Ap{oo) is 
completely invariant. 

We call z a fixed point of p if p(z) = z and a periodic point if p"(z) = z, 
for some n G N, in which case the finite sequence (zq, zi, ..., Zn_i), where zq = z, 
zi = p(- 2 o)) is the cycle of the periodic point z. We say that z is an attract¬ 
ing fixed point if |p^(-z)| < 1, a repelling fixed point if |p^(z:)| > 1, and a neutral 
fixed point if |p^(^)| = 1. Generalising, we say that a periodic point z is attract¬ 
ing/repelling/neutral if |P'(z)| < l/> l/= 1, where P = p". By the chain rule, 
P'(z) = p'(zo)...p'(z„_i), where zq = z and Zj := p(zj_i), j = l,...,n- 1. 

The value 7 = Pf z) is the multiplier of the neutral periodic point z (clearly 
I 7 I = 1). If z is a neutral periodic point with multiplier 7 we say that it is rationally 
neutral if 7 ^ = 1 for some iV G N, otherwise we say that it is irrationally neutral. 
We call z a critical point of p if p'{z) = 0. 
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If w is an attracting periodic point we denote by Ap{w) the hasin of attraction 
of the cycle C = {zq, Zn-i} of z, by which we mean Ap{w) := {z G C : 
d{p'^{z),C) —)■ 0 as n —>■ oo}. Here, for S' C C and z gC, d{z,S) := inf^jg^ \z — w\. 
It is easy to see that Ap{w) contains some neighbourhood of C, and hence that 
Ap(w) is open. 

We will make use of standard properties of the Julia set captured in the 
following theorem. We recall that a family F of analytic functions is normal at 
a point z G C if, in some fixed neighbourhood N of z, each f G F is analytic 
and every sequence drawn from F has a subsequence that is convergent uniformly 
either to some analytic function or to oo. 


Theorem 2.7. [Ill Summary 14.12] J{p) is compact with no isolated points, is 
uncountable, and has empty interior. J{p) is completely invariant, J{p) = J(jF), 
for every n G N, 

J{p) = {z G C : the family {p^,p '^,...} is not normal at z}, (2-9) 


J(p) is the closure of the repelling periodic points of p, and, for all except at most 
one z G C, 

OO 

•^(P) C IJ p-"({z}). (2.10) 

n—1 


The Fatou set F{p) has one unbounded component U. It follows from (2.9) 
that U C Hp(oo); indeed, U = Hp(oo) as a consequence of the maximum principle 
[2]. It may happen that this is the only component of F{p) so that Ap{oo) = F{p). 
This is the case if fc = (1,1) andp{z) = pk{z) = z'^ — 2, when iG(p) = J{p) = [—2, 2] 
[21 p. 55]. If F{p) has more than one component it either has two components (for 
example, if /c = (—1,1) and p{z) = Pk{z) = z^, when K{p) = D and J{p) = 9D), 
or infinitely many components m Theorem IV. 1.2]. It follows from ( |2.10| ) that 
J{p) = dAp(oo) and J{p) = dAp(w) if w is an attracting fixed point or periodic 
point [21 Theorem III.2.1]. Arguing similarly [21 Theorem 1.7], J{p) = dFgip), 
where Fb{p) '■= int{K{p)) = F{p) \ Ap{oo), so that K{p) = Fb{p). 

Because J(p) is completely invariant and p is an open map, the image V = 
p(U) of any component C/ of F(p) is also a component of F(p). Now consider 
the orbit of U, i.e. {p{U))^=i. The following statement of possible behaviours is 
essentially Sullivan’s theorem [21 pp. 69-71]. 


Theorem 2.8. Let U be a component of F{p). Then one of the following cases 
holds: 

i) p^{U) = U, for some n Gf^, in which case we call U a periodic component 
of F{p), and call the smallest n for which p^{U) = U the period of U. (If n = 1, 
when U is invariant, we also term U a fixed component of F(p).) 

a) p'^{U) is a periodic component of F(p) for some r G N, in which case we 
say that U is a preperiodic component of F(p). 


The above theorem makes clear that the orbit of every component U of F(p) 
enters a periodic cycle after a finite number of steps. To understand the eventual 
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fate under iterations of p of the components of the Fatou set it is helpful to 
understand the possible behaviours of a periodic component. This is achieved in 
the classification theorem (e.g. 0)- To state this theorem we introduce further 
terminology. Let us call a fixed component U of F{p) a parabolic component if 
there exists a neutral fixed point w € dU with multiplier 1 such that the orbit 
of every z G U converges to w. Call a fixed component U of F{p) a Siegel disk 
if p is conjugate to an irrational rotation on U, which means that there exists a 
conformal mapping ip : U ^ V with 0 G V and an irrational 0 S K such that 

(p{p{z)) = g{p{z)) =-fp{z), zGU, (2.11) 

where g{w) = 'yw, w G V, and 7 = exp(27ri0). It is easy to see that, for w G U, 
p(w) = w iS w = and p'{w) = 7. Thus every Siegel disk contains a unique 

irrationally neutral fixed point (the Siegel disk fixed point). 

Theorem 2.9. Classification Theorem [21 Theorem IV.2.1]. If U is a periodic com¬ 
ponent of F(jp) with period n S N, m which case U is also a component of 
F{p^) = F{p), then exactly one of the following 

a) U contains an attracting periodic point w which is an attracting fixed point 
of p"^, and U C Ap{w); 

b) U is a parabolic component of F{p'^); 

c) U is a Siegel disk component of FijF). 

The following proposition relates the above cases to critical points of p (see 
[21 Theorems III 2.2 and 2.3, pp. 83-84]): 

Proposition 2.10. IfU is a periodic component of F{p) with period n then either: 
(i) U is a parabolic component of F{p'^) or contains an attracting periodic point, 
in which case U'^^ip'^fU) contains a critical point of p; or (ii) U is a Siegel disk 
component of F{p'^) and there is a critical point w of p such that the orbit of w is 
dense in dU. 


The following proposition will do the work for us in Section [^ 


Proposition 2.11. Suppose that S' C C is bounded, open and simply-connected, that 
T C S is closed, and that the orbit of every critical point in K(p) eventually lies 
in T. Then _ 

00 

K{p) CG:=\J p-{S). 

n—1 


Proof. That z e G if z G J{p) follows from ( |2.10 ). li z G Fb(jP) then, by Theo¬ 
rems |2.8| and [2.91 after a finite number of iterations the orbit of z is in a periodic 
component of F{p) that is parabolic, part of the domain of attraction of an at¬ 
tracting periodic point, or is a Siegel disk. In the first two cases it follows that 
d(p”(z),C) —)■ 0 as n —)• 00 for some cycle G, but also d{p'^{w),G) -G 0 for some 
critical point w by Proposition 2.10[ This last implies that C C T is non-empty. 


^The result as stated for rational p in [2] gives a 4th option, that t/ is a Herman ring component 
of This option is excluded if p is a polynomial |21l p. 166]. 
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and so p^{z) G S for some n. In the case that the orbit of z is eventually in a 
Siegel disk then also p"(z) G S for some n for, if the orbit of every critical point 
w G K(p) is eventually in T, it follows that the boundary of every Siegel disk is in 
T, and (as S is simply connected) that every Siegel disk is in S. □ 


Previous upper bounds on E. We have noted above that, if c S 11 is pseudo- 
ergodic, then E = specie C W{Ac) = A, given by (1.31. Similarly, the spectrum 
of is contained in the closure of its numerical range, so tha 10 

E C {±'/z : z G spec (A^)} C N 2 := {±-\/2 : z G W{A^)}. (2.12) 


Hagger [15] introduces a new, general method for computing numerical ranges of 
inhnite tridiagonal matrices via the Schur test, which he applies to computing the 
numerical range of A^ (by expressing it as the direct sum of tridiagonal matrices). 
These calculations show that E C IV 2 C A; indeed, the calculations in m imply 
that N 2 = {rexp(i0) : 0 < r < p{0), 0 < 0 < 27r}, where p G C'(M) is even and 
periodic with period 7 r/ 2 , given explicitly on [— 7 r/ 4 , 7 r/ 4 ] by 


P(0) 


\/2, 7r/6 < |6*| < 7r/4, 

2/-y/cos2d + \/3 I sin20|, |0| < 7r/6. 


(2.13) 


(2.14) 


By comparison, in polar form, 

W{Ac) = A = {re*® : 0 < r < 2/(| cos0| + | sin0|), 0 < 0 < 27r}. 

Figure [^includes a visualisation of A and A 2 . 

The bound \2A2 ] , expressed concretely through (2.13), is the sharpest explicit 
upper bound on E obtained to date. It implies that E is not convex since we also 
know (see (3.1) below) that ±2, ±2i, and ±1 ± i are all in E. 

A different family of upper bounds was established in [3] (and see 0) , ex¬ 
pressed in terms of pseudospectra. For a square matrix A of order n and e > 
0 let speCgA denote the e-pseudospectrum of A (with respect to the 2-norm), 
i.e., speCgA := {A G C : ||(A — A /„)“^||2 > e“^}, with the understanding that 
II (A — A /)“^||2 = 00 if a is an eigenvalue, so that spec A c spec^A. (Here || • II 2 is 
the operator norm of a linear mapping on C" equipped with the 2-norm.) Analo¬ 
gously to (1.9), for e > 0 and n G N, let 

o-n.£ := IJ spec.A, (2.15) 


which is the union of the pseudospectra of 2"“^ distinct matrices. Then it is shown 
in [ 3 ] that 

E* := \ E as n — 00 , (2-16) 

where := 4sin0„ < 27r/(n -I- 2) and 0„ is the unique solution in (7r/(2n -|- 
6 ), tt/ (2n + 4)] of the equation 2 cos((n -|- 1)0) = cos((n — 1)0). 


^Equation ( |2.12[ l is the idea behind higher order numerical ranges; indeed, where p is the poly¬ 
nomial p(A) = A^, N 2 is Num(p, Ac) in the notation of CqI p. 278], so that N 2 is a superset of 
the second order numerical range. 
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Clearly, {E* : n e N} is a convergent family of upper bounds for E that is 
in principle computable; deciding whether A G E* requires only computation of 
smallest singular values of n x n matrices (see 0 (39)]). Explicitly E]) = 2D, and 
E* is plotted for n = 6,12,18 in [4j. But for these values E* D A, and computing 
E* for larger n is challenging, requiring computation of the smallest singular value 
of 2”“^ matrices of order n to decide whether a particular A G E*. Substantial 
numerical calculations in [3] established that 1.5 + 0.5i ^ E 34 , providing the first 
proof that E is a strict subset of A, this conhrmed now by the simple explicit 


bound (2.12) and (2.13). 


3. Lower Bounds on S and Symmetries of S and Stt 

Complementing the upper bounds on E that we have just discussed, lower bounds 


on E have been obtained by two methods of argument. The first is that (1.2) tells 


us that specHf, C E for every 6 G fl. In particular this holds in the case when b is 


periodic, when the spectrum of Ab is given explicitly by Lemmas 2.2 and |2.51 so 
that, as observed in the introduction, 

:= y specHP®''cE. 

fcG{±l}" 

Explicitly [U Lemma 2.6], in particular, 

TTi = [—2, 2] U i[—2, 2] and 7 r 2 = tti U {x ± ix : —1 < x < 1}. (3.1) 

In the introduction we have defined tToo := and have termed E^^. := W^, 

also a subset of E since E is closed, the periodic part of E. We have also recalled 
the conjecture of [3] that E^r = E. Let 


n. 


-u 


C TToo C E^ C E. 


Then it follows from Lemma 12.11 that 


n„ 


E^ as n —)■ 00 . (3.2) 

If, as conjectured, Ejr = E, then (3.2) complements ( 2.16[ ); together they sandwich 
E by convergent sequences of upper (E*) and lower (n„) bounds that can both 
be computed by calculating eigenvalues of n x n matrices. Figures and [^include 
visualisations of ttso, indistinguishable by eye from Hag, but note that the solid 
appearance of TTgo, which is the union of a large but finite number of analytic arcs, 
is illusory. See 00 for visualisations of 7 r„ for a range of n, suggestive that the 
convergence ( |3.2[ ) is approximately achieved by n = 30. 

The same method of argument to obtain lower bounds was used in 0 , 
where a special sequence b G was constructed with the property that spec Ab D 


D, so that, by (1.2), D C E. The stronger result (1.5), that this new lower bound 
on E is in fact also a subset of E.^., was shown in 0, via a second method of 
argument for constructing lower bounds, based on surprising symmetries of E and 
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Stt- We will spell out in a moment these symmetries (one of these described first 
in [5] , the whole infinite family in m), which will be both a main subject of study 
and a main tool for argument in this paper. But first we note more straightforward 
but important symmetries. In this lemma and throughout A denotes the complex 
conjugate of A € C. 

Lemma 3.1. [H Lemma 3.4] (and see [TH], 0 Lemma 4\). All of Tin, an, S-tt, and 
S are invariant with respect to the maps A i—>■ iA and A i—>■ A, and so are invariant 
under the dihedral symmetry group D 2 generated by these two maps. 


To expand on the brief discussion in the introduction, m proves the exis¬ 
tence of an infinite set S of monic polynomials of degree > 2, this set defined 
constructively in the following theorem, such that the elements p G S are symme¬ 
tries of TToo and S in the sense that (3.3) below holds. 


Theorem 3.2. [T7] Let S denote the set of those polynomials pk, defined by ( |2.5[ ), 
with k = (fci, ...,kn) G {±1}" for some n>2, for which it holds that: (i) fc„_i = 
— 1 and kn = 1; (ii) Pk = P%, where k G {il}" is the vector identical to k but with 
the last two entries interchanged, so that kn-i = 1 and kn = —1. Then 

T, C p{T,) and p~^{noo) C TToo, (3.3) 

for all p G S. 


We will call S Hagger’s set of polynomial symmetries for S. 

We remark that ii p G S then it follows from (3.3), by taking closures and 
recalling that p is continuous, that also 


p ^(S,r) C and p ^(int(S,r)) C int(I]^). 


(3.4) 


We note also that p ^(tToo) C tToo implies that tToo C p{tToo), but not vice versa, 
and that S C p{T,) iff 

for all Age. 


Further, we note that it was shown earlier in that (3.3) holds for the particular 
case p(A) = A^ (this the only element of S of degree 2, see Table JB; in [S] it was 
also shown, as an immediate consequence of (3.3) and Lemma 3.1 that 

p-i(E) CE, 


for p(A) = A^. Whether this last inclusion holds in fact for all p G 5 is an open 
problem. 

Our first result is a much more explicit characterisation of S. 


Proposition 3.3. The set S is given by S = {pk '■ k G K-}, where 1C consists of those 
vectors k = (ki, ..., kn) G {±1}" with n >2, for which: (i) kn-i = —1 and kn = 1; 
and (ii) n = 2, or n> 3 and kj = kn-j-i, for 1 < j < n — 2, so that {ki, ..., kn- 2 ) 
is a palindrome. Moreover, if k G 1C, then 

Pk{\) = A(3'fc(i:„_2)(A). 


(3.5) 
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Proof. It is clear from Theorem 3.2 that what we have to prove is that, iik G {±1}" 
with n > 2 and kn-i = —1, kn = 1 , then p/j = pj if n = 2 or 3; further, if n > 4, 
then pk = pj iff (/ci,..., k„_ 2 ) is a palindrome. 

li k G {±1}" with n>2 and kn-i = — 1 , = 1, then, from (2.6) and (2.3), 

Pk{^) ^nQk{2:n—2)i^) 

= Mk(l:n-2)W + Qk(l:n-3)W ~ <lk{2:n-2)W ■ 

Thus, if n = 2 or 3, or n > 4 and (ki, ...,kn-2) is a palindrome, PkW = 
Mk(i-n-2)W since qk{i-.n-3)W = ' 7 fc( 2 :n- 2 )(A), this a consequence of the defini¬ 
tions (2.2) in the cases n = 2 and 3, of Lemma 2.4 in the case n > 4. Similarly, 
Pfc(A) = Agfe(i,„_ 2 )(A), so that pk = p%. 

Conversely, assume that k G {±1}" with n > 4, kn-i = — 1 , kn = 1, and 
Pk = Pk- To show that (ki ,..., kn-2) is a palindrome we need to show that kj = 
kn-j-i, for 1 < j < (n — 2)/2. Using (2.6) and then ( |2.3| ), we see that 

0 =Pfc(A) -pg(A) = qk(l-.n-l)W - qk(i,n_i)W -‘2qki2-.n-2)W 

= 2gfc(i.„_3)(A) — 2qfc(2:„_2)(A). 

Thus qk{l:n—3) ~ qk{2-.n—2)- Hut, if qk(^j:n—j — 2) ~ qk{j + l:n—j — l) S-ud 1 ^ J ^ {ri 
2)/2, then, applying ([^, 


0 = 


j —2) (A) qk(^j.\-3:n—j—l){^') 

A^7fc(j-t-l:n—j —2) (A) kjqk(^j+2:n—j — 2)i^) 

(A^7fc(j;-t-l:n—j — 2 ) (A) kn—j—iqk(j + l:n—j—3){^)) 

^j^fc(j-t- 2 :n—j — 2 ) (A) “1“ kn—j — iqk(j + l:n—j — 3){^'}- 


As this holds for all A and, by Lemma |2.3| and (2.2), qk(j+ 2 :n-j- 2 ) and qk(j+i-.n-j- 3 ) 
are both monic polynomials of degree n — 2j — 2, it follows first that kj = kn-j-i 
and then that qk{j+i,n-]-3) = qk(j+2-.n-i-2)- Thus that kj = kn-j-i for 1 < j < 
{n — 2)/2 follows by induction on j. □ 


The following corollary is immediate from (3.5) and Lemma 2.3 


Corollary 3.4. Suppose that n > 2, k G {±1}", and pk G S. Then, as A 
PkW = ±A -I- O(A^) ifn is odd, while PkW = O(A^) if n is even. 


0 , 


Let us denote by Pm the polynomialp^ when k has length m > 2, km-i = —1, 
km = 1, and all other entries are I’s. It is convenient also to define Ti(A) = A. 
Clearly, as a consequence of the above proposition. Pm G S for m> 2 (that these 
particular polynomials are in S was observed earlier in El)- We will write down 
shortly an explicit formula for Pm in terms of Chebyshev polynomials of the 2nd 
kind. Recall that Un{x), the Chebychev polynomial of the 2nd kind of degree n, 
is defined by [T] Uo{x) := I, Ui{x) := 2x, and Un+i{x) := 2xUn{x) — Un-i{x), for 
n e N. 


Lemma 3.5. For m gN, PmW = A[/m-i(A/2). 


Proof. This follows easily by induction from (3.5) and (2.3). 


□ 
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We note that, using the standard trigonometric representations for the Cheby- 
chev polynomials [1], for to S N, 

Pm(2 cos 0) = 2 cos 9Um-i{cos 0) = 2cot0sinTO0 =: rm{0)- (3.6) 

A similar representation in terms of hyperbolic functions can be given for the 
polynomial pfc when k has length 2 to— 1 and kj = (—1)-^; we denote this polynomial 
by Qm- Clearly, for to > 2, € 5 by Proposition [3]^ and Qm is an odd function 

by Lemma |2.5| The proof of the following lemma, like that of Lemma |3.5[ is a 
straightforward induction that we leave to the reader. 


Lemma 3.6. Qi(A) = A, Q 2 (A) = A^ + A, and Qm+i(A) = + Qm-i(A), 

for TO > 2. Moreover, for to G N and 9 > Q, 


Qm (^V2 sinhd^ 


y/ 2 sinh 9 
V 2 sinh 9 


sinh(TO0) + cosh((TO — l)d) 

cosh 9 ’ 

cosh(TO0) + sinh((TO — l)d) 

cosh 9 ’ 


if m is even, 
if TO is odd. 


The following lemma leads, using Lemmas 3.5 and |3.6[ to explicit formulae 
for other polynomials in S. For example, if denotes the polynomial pk when k 
has length to > 2, km-i = ~1, km = 1, and all other entries are —I’s, then, by 
Lemmas 13.51 and 13.71 


P,;(A) = i-^PmiiX) = i^-^XUm-iiiX/2). 


(3.7) 


Lemma 3.7. If k G {±1}" and pj. G S, then p-k & S and p-k{X) = i "pfc(iA). 


Proof. Suppose that k G {±1}" and pk G S. li n = 2, then k = —k and p-k = 
p^ = Pk G S. li n > S, defining £ G {±1}" by £n-i = —1, in = 1, and ij = —kj, 
for j = 1,..., n — 2, pg G S hy Propositi on |3.3[ so that p-k = pf = pt G S. That 
P-k{X) = i~"'Pfe(iA) comes from Lemma [zep □ 


3.3 


implies that there are precisely 2 




vectors 


We note that Proposition 
of length n in /C, so that there are between 1 and 2\^'^~^ polynomials of degree 
n in 5, as conjectured in m- Note, however, that there may be more than one 
k G 1 C that induce the same polynomial pk G S. For example, Pfc(A) = A® — A^ for 
k = (—1,1,1,-1, —1,1), and, defining £ = (1, —1, —1,1, —1,1) and using Lemma 
also pe{X) = Pi{X) = p-k{X) = —pfc(iA) = A® — A^. In Table (cf. [17]) we 


3.7 


tabulate all the polynomials in S of degree < 6. 

If p, <7 G S, so that p and q are polynomial symmetries of S in the sense 


that (3.3) holds, then also the composition p o g is a polynomial symmetry of S in 
the same sense. But note from Table that, while P3 o P2 G S, none of P2 o P2, 
P2 o P2 o P2, Q 20 P 2 , P 2 ° P 3 , or P2 o Q2 are in S. Thus S does not contain all 
polynomial symmetries of E, but whether there are polynomial symmetries that 
are not either in S or else compositions of elements of S is an open question. 

We finish this section by showing in subsection |3.1| the surprising result that 
S is large enough that we can reconstruct the whole of from the polynomials 
Pk G S. This result will in turn be key to the proof of our main theorem in Section 
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k 

PkW 

(-1,1) 

= P2(A) 

(1,-1,1) 

A^ - A = P3(A) 

(-1,-1,!) 

A^ + A = Q2(A) = P3*(A) 

(1,1,-1,1) 

A^ - 2A^ = P4(A) 

(-1,-1,-1,1) 

A4 + 2A^ =P|(A) 

(1,1,1,-1,1) 

A’^ - 3A^ + A = P5(A) 

(1,-1,1,-1,1) 

A^-A^ + A = -iQ3(iA) 

(-1,1,-1,-1,1) 

A^ + A^ + A = Q3(A) 

(-1,-1,-1,-1,1) 

A'> + 3A^ + A = P5*(A) 

(1,1,1,1,-1,1) 

A'’-4A4 + 3A" =P6(A) 

(1,-1,-1,1,-1,1) 

A«-A" =P3(P2(A)) 

(-1,-1,-1,-1,-1,1) 

A‘’+4A4 + 3A^ = P*(A) 


Table 1. The elements pk G S oi degree < 6. 


1^ Then in subsection 3.2 we use that (3.3) holds for the polynomials in Table 
to obtain new explicit lower bounds on including that I.ID C C E. 


3.1. Connecting eigenvalues of finite matrices and polynomial symmetries of E 

Recall from Proposition |3.3| that S = {pk ■ k S /C}, let K := U„gN{±l}", and 
define 

■= U spec HP®" c TToo = U specHf". (3.8) 

fce/c k^K 

The following result seems rather surprising, given that K, is much smaller than K 
in the sense that there are precisely vectors of length n in /C but 2” in K. 

Theorem 3.8. CToo C so that is dense in E,r and := tToo = CToo = 

Proof. We will show in Proposition |3.9| below that, for n > 2, 

CTn C 7rf„_^2 := U specHP®" C 7r2„+2, 

^^^2ti.+2 

where, for m > 2, Km denotes the set of those vectors in K that have length m. 
Since also CTi = {0} C tt^, for every m G N (e.g. [H Lemma 2.10]), this implies 
that (Too C TT^, which implies that tt^ is dense in E,r since (Too is dense in tToo [IHl 
Theorem 1 ], and the result follows. □ 


The key step in the proof of the above theorem is the following refinement of 
Theorem 4.1 in [4], which uses our new characterisation, Proposition 3.3 of S. 


Proposition 3.9. Suppose a, b,c,dG {±1} and k G {±1}”, for some n>2, and let 
k := (fci,..., fc„_i). Then 

spec C spec HP®", for £ = (fc', a, 6, fc, c, d) G {il}^"''"^, (3.9) 

where k' = kJn-i = (kn-i, ■■■,ki). Further, specH^"^ C 7rf„+2- 
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Proof. The proof modifies [H Theorem 4.1] where the same result is proved for 
the special case that a = c = —l,5 = d= l. Following that proof, suppose that 
A is an eigenvalue of with corresponding eigenvector x, let x := J„x and 


R") ■— J A 


(n) 


Jji 


SO that A^r'^x = Ax, and set 


B := 



All 


V 


-1 ■■■ / 


where [oj marks the entry at position (—1, —1). Then B is & bi-infinite tridiagonal 
matrix with zeros on the main diagonal, and with each of the first sub- and super¬ 
diagonals a vector in id that is periodic with period 2n + 2. Define x € the 
space of bounded, complex-valued sequences : Z —>■ C, by 

X := (...,0,x^,[^, x^,0,x^,0,x^, ...)^, 

where marks the entry x_i. Then it is easy to see that Bx = Ax, so that 
A G spec bQ Further, by a simple gauge transformation [H Lemma 3.2], specB = 
spec where £ is given by (3.9). 


We have shown that spec C spec But, choosing in particular a = h 
and c = —l,d=l,we see from Proposition ) 
spec^P” C 7r|„+2- 


3.3 


that ^ G lC 2 n+ 2 i so that spec c 

□ 


Remark 3.10. Proposition 3.9 implies that spec^l"^ C specA?®“^ for 16 different 


vectors £ G {±1}", corresponding to different choices of a,b,c,d G {±1}- Some of 


these vectors correspond to the same polynomial pg and hence, by (2.8), to the same 


spectrum specIn particular, if a = 6 = ±1, then the choices c = —d = 1 and 
c = —d = — 1 lead to the same polynomial by Proposition |3.3| and the definition 
of S. But, if a ^ again by Proposition |3.3| and the definition of S, the choices 
c = —d = 1 and c = —d = —1 must lead to different polynomials, and neither 
of these polynomials can be in S. On the other hand, as observed already in 
the proof of Proposition |3.9[ the choices a = b and c = —d lead to an i G K, 


“^Clearly A is in the spectrum of B as an operator on f°°, but this is the same as the t^-spectrum 
from general results on band operators, e.g. |20]. 












Feinberg-Zee Random Hopping Matrix 


17 



Figure 1. An illustration of Proposition |3.9| in the case n = 3, 
ki = k 2 = 1. The red circles indicate the eigenvalues, 0 and ±-\/2, 
of The black lines are the spectra of for the different 
choices of £ defined by (3.91. In this case there are 7 distinct 
polynomials pi and 7 associated distinct spectra specA^®'^, each 


(S') 

of which contains the eigenvalues of A^ b One cannot see all the 
spectra as separate curves because some of them overlap. 


and so to a polynomial pg G S. Thus Proposition |3.9| implies that there are at 
least three distinct polynomials pi such that specA^"^ C spec A^®"^. Figure plots 


specA^,"^ and specAP®"^, for the vectors defined by (3.9), in the case that n = 3 


and ki = k 2 = I- For other plots of the spectra of finite and periodic Feinberg-Zee 
matrices, and the interrelation of these spectra, see [HElIl]. 


3.2. Explicit lower bounds for Ett 

As noted in that the polynomials p G S satisfy ( |3.3[ ) gives us a tool to compute 
explicit lower bounds on E.^.. Indeed m shows visualisations of p "(D), for several 
p G S and n G N, and visualisations of unions of p“"(D) for p varying over some 
finite subset of S. Since, by (1.5), D C E^, it follows from (3.3) that all these sets 
are subsets of E,r- 

The study in m contains visualisations of subsets of E^ as just described, 
but no associated analytical calculations. Complementing the study in we 
make explicit calculations in this section that illustrate the use of the polynomial 
symmetries to compute explicit formulae for regions of the complex plane that are 
subsets of Ett, adding to the already known fact (1.5) that D C E^^.. 
















18 


Simon Chandler-Wilde and Raffael Hagger 


Our first lemma and corollary give explicit formulae for p“^(D) when p(A) = 
PsW = — X and when p{X) = Q 2 {X) = + A. These formulae are expressed 

in terms of the function s G C°°[—1,1], where, for — 1 < t < 1, s(t) > 0 is defined 
as the largest positive solution of 

/(s) := — 2ts^ + s = 1. 


If —1 < t < 1/2, this is the unique solution in (0,1) (on which interval f'(s) > 0), 
while if 1/2 < t < 1 it is the unique solution in [1, 2) (on which interval f'(s) > 0). 
Explicitly, for — 1 < t < 1, 


s'(t) 


2-s(t) + (s(t))3 


> 0 , 


and s(l/2) = 1, s(l) « 1.75488. 


Lemma 3.11. If p{X) := P 3 (A) = A^ — A, then 

p-i(D) = E-.= {re‘® : 0 < r < 5'(6»),0 < 0 < 2tt}, 

where S{9) := a/s (cos2d), for d € M. S' € C°°(K) and is even and periodic with 
period n. In the interval [0,7r] the only stationary points of S are global maxima 
at 0 and tt, with S(0) = S(7r) = -/s/T) « 1.32472, and a global minimum at 'k 12. 
Further, S is strictly decreasing on [0,7r/2] and S{9) >1 m [0,7r] iff 0 < 0 < tt/G 
or 57r/6 < d < tt, with equality iff 9 = tt/G or Stt/G. 


Proof. If A = re‘® (with r > 0, d S K), then |p(A)p = jA^ — Ap = jr^e^'® — rp = 
r® —2r"‘cos(2d) + r^. It is straightforward to show that |p(A)| < 1 iff 0 < r < S(d). 
The properties of S claimed follow easily from the properties of s stated above. □ 


Corollary 3.12. If p{X) := Q 2 (X) = A^ + A, then 

p-i(D) =iE = {re'^ : 0 < r < S(9 - 7r/2), 0 < d < 27r}. 


Proof. This is clear from Lemma 3.11 and the observation that Q 2 (A) = iP 3 (iA). 

□ 


Since P'i,Q 2 G S, it follows from the above lemma and corollary and (3.4) 
that EUiE C int(E,r)- But this implies by ( |3.4[ ), since T 2 (A) = A^ is also in S, that 
also ViE C int(S,r), where, for S C C, VS := {X G C : X'^ G S}. In particular, 

Wi := {re'® : 0 < r < S{9),—-k/Q < 9 < tt/G} C E C int(E,r) and (3.10) 
IT 2 := jre'® : 0 < r < a / S(2d — 7r/2), tt/G < d < 7r/3} C ViE C int(E,r)- (3-11) 


It is easy to check that 
3 

y i’"(ITi U IT 2 ) = E U iE U /iE C int(E^). (3.12) 

m—0 

Next note from Table that p G S where p{X) := A® — A^ + A factorises as 
p(A) = A(A - e'’^/®)(A + e'’^'''®)(A - e-'’^/®)(A + e'^’^/®). 
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Figure 2. A plot showing Wi (green), W 2 (red), + tjD 

(blue), and their rotations by multiples of 7r/2. The union of the 
green, red, and blue regions is VF C defined by (3.14[). W 


contains I.ID, indicated by the black circle, see Proposition |3. 13 
In the background in grey one can see tt^q C The dotted and 
dashed-dotted curves are the boundaries of A and A 2 , respec¬ 
tively, defined by (1.3) and (2.12), with A D A 2 D E D E^. 


Thus, for A = exp(±i7r/6) -I- w with |ty| < e, 

b(A)| < (l-l-e)e(2-|-e)(2sin(7r/6)-|-e)(2cos(7r/6)+e) = e(l-|-e)^(v^-|-e)(2-|-e) =: g(e). 

Let ry « 0.174744 be the unique positive solution of g(e) = 1. Clearly |p(A)| < 1 if 
A = exp(±i7r/6) -I- w, with |r(;| < ry, so that 

exp(±i7r/6) + lyD C p“^(D) C int(E^). (3.13) 

We have shown most of the following proposition that extends to a region 
W D I.ID (illustrated in Figure]^ the part of the complex plane that is known 
to consist of interior points of E^, making explicit implications of the polynomial 
symmetries 5 of E. Before m and the current paper the most that was known 
explicitly was that D C int(E,r)- 
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Proposition 3.13. 

3 

l.m CW:= y i'" (Wi UW 2 U + t]B) U 

m—0 


C int(E^). 

(3.14) 


Proof. That W C int(E^) is (3.12) combined with (3.13) and Lemma 3.1 It is 
easy to see that W is invariant with respect to the maps A 1 —> iA and A 1 —A, 
and so is invariant under the dihedral symmetry group D 2 generated by these two 
maps. Thus to complete the proof it is sufficient to show that z := re'® G W for 
0 < r < 1.1 and 0 < 0 < 7r/4. 


Now since, by Lemma 3.11 S(0) > 1 for —ttI& <d<7r/6, zG Wi U W 2 

S is even and 


for 0 < r < 1, 0 < 0 < 7r/4. Further, since, by Lemma 3.11 


is increasing on [—7r/6, 0], S{9) > S{Tr/8) = J s(-\/2/2), for 0 < 0 < tt/S, and 


a/ S{29 — 7r/2) > -y/5'(7r/12) = (s(-\/3/2)) ^ , for 57r/24 < 9 < 7r/4. But -s/S > 
1.73, so that /(1.5) < 0.9825 for t = which implies that s(\/3/2) > 1.5, 

so that (s(-\/3/2))^^^ > 1.1. Similarly, •\/2/2 > 0.705 and s(0.705) = 1.25, so that 
J s(-\/2/2) > 1/1-25 > 1.1. Thus z G Wi for 0 < r < 1.1 and 0 < 0 < tt/S, while 


z G IT 2 for 0 < r < 1.1 and 57r/24 < 9 < -k j4:. 

To conclude that I.ID C W, it remains to show that z G W ior 1 < r < 1.1 
and |7r/6 — 9\ < 7r/24. But it is easy to check that, for these ranges of r and 9, 
z G exp(i7r/6) + ryD provided cos(7r/24) + ^cos^(7r/24) + — 1 > 1.1. But this 

last inequality holds since cos(7r/24) = f y^2 + -\/2 + > 0.991 and g(0.174) < 1 

so that rj > 0.174. □ 


3.13 


extended to a region W D I.ID the part of the 


4. Interior points of 

We have just, in Proposition 
complex plane that is known to consist of interior points of In this section we 
explore the relationship between and its interior further. We show first of all, 
using (3.4) and that P„ G 5 for every n > 2, that [0,2) C int(S^). Next we use 
this result to show that, for every n > 2, all but finitely many points in 7rf are 
interior points of E^^.. Finally, we prove, using Theorem |3.8[ that E^ is the closure 
of its interior. If indeed it can be shown, as conjectured in [3], that E^ = E, then 
the result will imply the truth of another conjecture in [3], that E is the closure 
of its interior. 

Our technique for establishing that [0,2) C int(E 7 r) will be to use that 
1,1)) C int(E,r), for every m > 3, this a particular instance of (3.4). 
This requires first a study of the real solutions of the equations Pm (A) = ±1 and 
their interlacing properties, which we now undertake. 
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From ( |3.6| ), 

Pm{2) = 2m, Pm(2cos( 7 r/m)) = 0, and Pm(2cos( 37 r/( 2 m))) = —2cot( 37 r/( 2 m)). 

This implies that the equation Pm{^) = 1 has a solution in (2 cos( 7 r/m), 2). Let 
A+ denote the largest solution in this interval. Further, if m > 5 then Pm (A) = — 1 
has a solution in (2 cos( 37 r/( 2 m)), 2) since —2 cot( 37 r/( 2 m)) < —1. For m > 4 let 
A“ denote the largest solution to Pm(A) = —1 in (0,2), which is in the interval 
( 2 cos( 37 r/( 2 m)), 2 ) if m > 5, while an explicit calculation gives that A 4 = 1. 

Throughout the following calculations we use the notation r„(0) from (3.6). 

Lemma 4.1. For m > A it holds that Pm is strictly increasing on (A“,2), that 
A“ < A+, and that —1 < Pm(A) < 1 for A“ < A < A+. 

Proof. Explicitly, P 4 (A) = XU‘i{X/2) = X'^ — 2A^, so that P 4 (A) = 4(A^ — 1) and 
these claims are clear for m = 4. 

Suppose now that m > 5. It follows by induction that, for n > 3, r„(0) is 
strictly decreasing on ( 0 , 7 r/n + 7 r/n^). For r^{9) = 2cos0(4cos^ 0 — 1) is strictly 
decreasing on ( 0 , 47 r/ 9 ) C ( 0 , 7 r/ 2 ) and, if this statement is true for some n > 3, 
then 

r„+i( 0 ) = 2 cot 0 sin((n + 1)0) = cos 0 (r„( 0 ) + 2 cosn 0 ) 
is strictly decreasing on ( 0 , 7 r/(n + 1) + 7 r/(n + 1)^) C ( 0 , 7 r/n). Further, 

rm(7r/m + 7r/m^) = —2cos(7r/m + 7r/m^) sin(7r/m)/sin(7r/m + 7r/m^) 


< - 


2mcos(7r/5 + 7r/25) 


< -10cos(67r/25)/6 < -1, 


m + 1 

since sina/sin6 > a/h for 0 < a < 6 < tt. As rm(0) = Pm(2cos0), these obser¬ 
vations imply that, on (2cos(7r/m -I- 'kI rnff 2), Pm is strictly increasing, and that 
A“ > 2cos(7r/m -|- 7r/m^). Thus A+ > A“ follows from the definitions of Xijf, and 
— 1 < Pm(A) < 1 for A“ < A < A+. □ 


As observed above it follows from (3.4) that P^^((—1,1)) C int(S^), for 
m > 2. Combining this observation with Lemma 4.1 and the fact that PsiX) = 
A3-As (-1,1) if-A+ < A< A+, we obtain the following corollary. 


Corollary 4.2. For m = 4, 5,..., (A,^ 


, A+) C int(S 7 r). Also, (—Ag ,A|^) C int(S^). 


By definition of A+, A+ —)■ 2 as m —)■ 00 . Thus the above corollary and the 
following lemma together imply that [0, 2) C int(E 7 r). 

Lemma 4.3. For m = 3,4,..., A “_|_4 < A+ < 


Proof Since P 3 (A) = A^ — A and P 4 (A) = A^ — 2A^, it is easy to see that 1 = A 4 < 
AJ < A 4 , so that the claimed result holds for m = 3. 

To see the result for m > 4 we will show, equivalently, that rf^_^_l > r+ > 
for m = 4,5,.... Here r+, for n e N, is the smallest solution of r„(0) = 1 
in ( 0 , 7 r/n), so that A+ = 2cosr+, while r~, for n = 4,5,..., denotes the smallest 
solution of r„(0) = —1 in ( 0 , 7 r/ 2 ), so that X~ = 2cosr“. 
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We have shown in the proof of Lemma 4.1 that, for m > 3, rm(P) is strictly 
decreasing on (0, Tr/m + Tr/m^), and that rj'm, ^ n jrr?) < —1, for m > 5, while 


for m > 5. 


rm(7r/TO) = 0 so that 

TT 

_ TT TT 


— 

2 ’ 


m 

m 

Similarly, for to > 2, 



rm(7r/"i - tt/to^) 

= 

2 cos(7r/m — 7rjm' 


> 

2 cos(7r/2 — 7r/4) = 

so that 

TT 

TT ^ vr 

— — <r <r — 

m 


m 


for TO > 2. 

These inequalities imply that, for to > 4, G (0,7r/TO + ttIvt?), and since 

is strictly decreasing on this interval and 

rm(r”+i) = 2cotr-+;^sin((TO +l)r-+i - r"+i) 

= cosr-_^i(-l-2cos((TO + l)r-+i)) <cosr-+i < 1, 

it follows that r+ < < 7r/(TO + 1) + 7r/(TO + 1)^. Since also r^+i is strictly 

decreasing on (0, tt/(to + 1) + TT/{m + 1)^) and 

rm+i{rt) = 2cotr+sin(TOr++r+) = cosr+(l + 2cos(TOr+)) < cosr+ < 1, 

since 7r/2 < TOr+ < tt, we see that also r+ > 

Corollary 4.4. (—2,2) Ui(—2, 2) C int(S^). 

and Lemma 

-i+i ^ 


□ 


Proof. From Corollary 4.2 and Lemma 4.3 it follows that [0, A^) C int(E,n.) and 
that [A+, C int(E^), for to > 3. Thus [0, A+) C int(S,r) for to > 3, and so 

[0,2) C int(I] 7 r) since A+ —?> 2 as to —>■ oo. Applying Lemma 3.1 we obtain the 
stated result. □ 


The following lemma follows immediately from Corollary 4.4 (2.8), and (3.4). 

Lemma 4.5. Suppose that k G so that k has length n and pk € S has degree n. 
Then all except at most 2n points in spec are interior points of Further, 
if X G spec then there exists a sequence (Am) C spec A^®'^ C int(E,n.) such that 
Am —t A as m ^ oo. 

As an example of the above lemma, suppose that k = (—1,1) S /C 2 . Then 
(see Table[^ Pfe(A) = A^ and, from (2.8), specAP®*^ = {cc ± ix : —1 < cc < 1}. 
There are precisely four points, ±l±i € spec A^®‘^\int(E^). These are not interior 
points of Sjr since they lie on the boundary of A D S D 

Combining the above lemma with Theorem |3.81 we obtain the last result of 
this section. 


Theorem 4.6. Ett is the closure of its interior. 
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Proof. Suppose A e Then, by Theorem 3.8 A is the limit of a sequence (A„) C 
TT^, and, by Lemma 4.5 for each n there exists S int(I] 7 r) such that |/in —A„| < 
n~^, so that —>■ A as n —)■ oo. □ 


5. Filled Julia sets in 

It was shown in m that, for every polynomial symmetry p € S, the corresponding 


Julia set J(p) satisfies J{p) C U{p) C S^r, where U{p) is defined by (1.7). (The 


argument in [IT] is that J{p) C U{p) by (2.10), and that U{p) C Stt by (3.4).) It 
was conjectured in m that also the filled Julia set K{p) C U{p) C S,r, for every 
p G iS. In this section we will first show by a counterexample that this conjecture 
is false; we will exhibit a p € 5 of degree 18 for which K{p) <f. U{p). However, we 
have no reason to doubt a modihed conjecture, that K{p) C S,r, for all p G S. 
And the main result of this section will be to prove that K{p) C for a large 
class of p S S, including p = Pm, for m > 2. 

Our first result is the claimed counterexample. 

Lemma 5.1. Let k = (1, —1,1,1,1, —1,1, —1, —1,1, —1,1,1,1, —1,1, —1,1), so that 
(by Proposition's^ pk G S, this polynomial given explicitly by 

Pfc(A) = A^® - 4A^® + 5Ai'‘ - 4A12 ^2^ 


Then K{pk) ^ C/(pfe). 

Proof. Let p = pfe. If we can find a /i ^ D that is an attracting fixed point of 
p, then, for all sufficiently small e > 0, := p + eD satisfies p{N) C N and 

TV n D = 0, so that N C K{p) and N r\U{p) = 0. Calculating in double-precision 
floating-point arithmetic in Matlab we see that A « 1.21544069 appears to be a 
fixed point of p, with 

p'(A) = 18Ai^ - 64Ai® -b 70Ai® - dSA^^ -b 70A^ - 64A^ -b 36A® - 16A® -b 2A « -0.69, 

so that this fixed point appears to be attracting. To put this on a rigorous footing 
we work in exact arithmetic to deduce, by the intermediate value theorem, that 
p(A) = A has a solution A* e (1.215,1.216), and that |p'(1.2155)| < 0.71. Then, 
noting that p"(A) = p+(A) — p_(A), where p+(A) = 306A^® -b 910A^^ -b 630A® -b 
180A^ -b 2 and P-(A) = OOOA^"^ -b 528A^° -b 448A® -b 48A^, we see that 

|p"(A)| < max{|p_(1.216) -p+(1.215)|, |p_(1.215) -p+(1.216)|} < 400 

for 1.215 < A < 1.216. But this implies that |p'(A*)| < |p'(1.2155)|-bO.0005x 400 < 
0.91, so that A* is an attracting fixed point. □ 

Numerical results suggest that amongst the polynomials p G 5 of degree < 20, 
there is only one other similar counterexample of a polynomial with an attracting 
fixed point outside the unit disk, the other example of degree 19. 

We turn now to positive results. Part of our argument will be to show, for 
every p G S, that {z : \z\ > 2} C AAoo), via the following lower bounds that 


follow immediately from Lemma |2.3[ (2.6) and (3.5). 
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Corollary 5.2. If k G {±1}", for some n G N, then |pfe(A)| > 2, for |A| > 2. If 
Pk G S, then |pfc(A)| > 2n, for |A| > 2. 

Corollary 5.3. Let p = pk, where k G {±1}", for some n G N. Then Ap{oo) D 
{z G C : |z| > 2}. If p G S, then Ap(oo) D {z G C : |z| > 2}. 


Proof Let z G C with |z| > 2. Then, by Corollary |5.2[ for some neighbourhood 
N of z, \pk{w)\ > 2 for w G N. Thus, and by Montel’s theorem [HI Theorem 
14.5], the family {p” : n G N} is normal at z. So z ^ J{p)y by (2.9). We have 
shown that J{p) n {z : jzj > 2} = 0, so that also K{p) n {z : jzj > 2} = 0 and 
^p(oo) = C \ iL(p) D {z : jzj > 2}. 

If p G 5 and |z| = 2 then, by Corollary 


5.2 


and so z G Ap(oo). Thus Ap(oo) D {z : jzj > 2]~ 


b(-z)| > 4 so that p{z) G Ap(oo) 

□ 


We remark that the bounds in Corollary |5.2| appear to be sharp. In particular, 
if fc = ( 1 , 1,..., 1) has length n > 2, we see from (2.5), ( |3.5[ ), and Lemma 3.5 that 
Pfc(2) = {7„{1) — 17„_2(1) = 2, since 17^(1) = m + 1 [I]. And we note that, if 
P = Prm for some m G N, then p{2) = Pm{2) = 2Um-i{I) = 2to. Finally, we recall 
that we have already noted that, for p = pk, with k = (1,1), i.e., p{z) = z^ — 2, 
the Julia set is J{p) = [—2, 2], so that Ap(oo) {z ; jzj > 2} for this p. 

The polynomial p{z) = z^ — 2 is an example where J{p) = K{p) so Fb{p) = 
K{p)\J{p) = %. The next lemma tells us that this does not happen, that K{p) is 
strictly larger than J(p), if p G S. 

Lemma 5.4. Fb(p) C U(p) is non-empty for p G S. 


Proof. If p G 5 is even then, by Lemma 2.5 and Corollary 3.4 p(0) = p'(0) = 0, 
so that 0 is an attracting fixed point. Clearly Ap(0) (which is non-empty) is a 
subset of U{p) n Fb{p). Similarly, by Lemma 2.5 and Corollary |3.4[ if p is odd 
then p(0) = 0 and p'(0) = ±1, so that 0 G J{p) is a rationally neutral fixed point 
and has a (non-empty) attracting region contained in Fb(j>) [H Section 11.5], this 
region clearly also in U{p). □ 


The above lemma and (3.4) imply that Fb{p) C D Fb{p) C U{p) is non¬ 


empty for all p G 5, in particular that Ap(0) C Fb{p) C 17(p) C if p is even. 
The main result of this section is the following criterion for the whole of Fb (p) to 
be contained in Ett. 

Theorem 5.5. Suppose that p G S , and that the critical points of p in K{p) have 
orbits that lie eventually in I.IDU (—2,2) Ui(—2, 2). Then K(p) C E.^.. 

Proof. Choose a and b with —2 < a < b < 2 such that K{p) DM. C [a,b] and 
K{p) n iM C i[a, &], this possible by Corollary 5.3 which says that the closed set 
K{p) C {z : \z\ < 2}. Set T = [a, b] U i[a, 6] U I.ID, and choose a simply-connected 
open set S such that T C S' C E^, this possible by Corollary |4.4| and Proposition 


3.13 By hypothesis, the orbits of the critical points in K{p) lie eventually in T. 
Thus the lemma follows from Proposition 2.11 and (3.4). □ 
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Figure 3 . Plot of 7130 C and (in blue) the fill ed Ju lia set K{p) 
in the casep(A) = Pl{\) = A^+2A^. By Corollary |5.6[ K{p) C 


As an example of application of this theorem, consider p G S given by (see 
Table p(A) = FKA) = A"^ + 2A ^. T his p has critical points 0 and ±i. Since 


p^(±i) = 3 it follows from Corollary 5.3 that ±i S Ap(oo), while 0 is a fixed point. 


Theorem |5.5| tells us that K{p), visualised in Figure]^ is contained in S^r. We note 
that, since all the critical points of p except the fixed point 0 are in Ap((X)), K(p) 


is not connected [U Theorem III 4.1] and, by Theorem |2.8| and Proposition 2.10 


Fb{p) = AL(0), which implies that K{p) C U{p). Further, recalling the discussion 
in Section!^ J{p) = dK{p) = 9Ap(0) = 9Ap(oo), and, since K{p) has more than 
one component, Fb(j>) has infinitely many components Theorem IV 1.2]. 

The above example is a particular instance of a more general result. It is 
straightforward to see that if p is a polynomial with zeros only on the real line, 
then all the critical points are also on the real line. Since, by Lemma 3.5 Pm (A) = 
A[/m-i(A/2), and all the zeros of the polynomial Um-i are real, it follows that all 
the zeros of Pm are real, so all its critical points are also real, and so the orbits of 


all the critical points are real. Further, by Corollary |5.3[ the orbits of the critical 
points in K{p) stay in (—2, 2). Likewise, as (see (3.7)) PmW = i”"*Pm(iA), all the 
critical points of P^ lie on iM, and so the orbits of these critical points are real if 
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m is even, pure imaginary if m is odd. Further, by Corollary |5.3| the orbits of the 


critical points in K{p) stay in (—2, 2) Ui(—2, 2). Applying Theorem 5.5 we obtain: 
Corollary 5.6. K{Pm) C Stt and K{P^) C Stt, for m> 2. 

Numerical experiments carried out for the polynomials p G S oi degree < 7 
(see Table and m Table 1]) appear to confirm that these polynomials satisfy 
the conditions of Theorem |5.5[ i.e., it appears for each polynomial p that the 
orbit of every critical point either diverges to infinity or is eventually in I.ID U 
(—2, 2) U i(— 2, 2). The same appears true for the polynomial p G 5 of degree 18 in 
Lemma 


5.1 


for which K{p) U{p). Thus it appears, from numerical evidence and 
Theorem 5.5 that K(p) C for these examples. These numerical experiments 
and Corollary |5.6| motivate a conjecture that K{p) C Stt for all p G S. 


6. Open Problems 

We finish this paper with a note of open problems regarding the spectrum of the 
Feinberg-Zee random hopping matrix, particularly problems that the above discus¬ 
sions have highlighted. We recall first that [3] made several conjectures regarding 
S. It was proved in |16j that but the following conjectures remain open: 

1. = E; 

2. E is the closure of its interior; 

3. E is simply connected; 

4. E has a fractal boundary. 

Of these conjectures, perhaps the first has the larger implications. Certainly, if E = 


Ett, then we have noted below (3.2) that we have constructed already convergent 


sequences of upper (E*) and lower (n„) bounds for E that can both be computed 
by calculating eigenvalues of n x n matrices. Further, if E = E^^., then the second 
of the above conjectures follows from Theorem |4.6[ 

The last three conjectures in the above list were prompted in large part by 
plots of 7r„ in |3], the plot of ttso reproduced in Figures]^ andIt is plausible that 
these plots, in view of (3.21, approximate Ejr. We see no clear route to establishing 


the third conjecture above. Regarding the fourth, we note that the existence of the 


set S of polynomial symmetries satisfying (3.3) suggests a self-similar structure to 


TToo and to E.^. and E and their boundaries. Further, m has shown that E^^. contains 
the Julia sets of all polynomials in S, and Proposition |5 . 5| and Corollary |5.6| show 
that E^ contains the filled Julia sets, many of which have fractal boundaries, of 
the polynomials in an infinite subset of S. 

Regarding these polynomial symmetries we make two further conjectures: 

5. K{p) C E,r for all p G iS; 

6. p“^(E) C E for all p G S. 

This last conjecture follows if E = E^r, by Theorem 3.2 from m- Further (see 
the discussion below (3.4)), it was shown in |S] that p“^(E) C E for the only 
polynomial of degree 2 in S, p(A) = A^. 
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The major subject of study and tool for argument in this paper has been 
Hagger’s set of polynomial symmetries S. We finish with one final open question 
raised immediately before Section [3T| 


7. Does S capture all the polynomial symmetries of S? Precisely, are there poly¬ 
nomial symmetries, satisfying (3.31, that are not either in S or compositions 
of elements of 5? 
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